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Abstract 

We simulate precision measurements of gauge-mediated supersymmetry breaking (GMSB) 
parameters at a 500 GeV e+e~ linear collider in the scenario where a neutralino is the next- 
to-lightest supersymmetric particle. Information on the supersymmetry breaking and the 
messenger sectors of the theory is extracted from the measured sparticle mass spectrum 
and neutralino lifetime. 
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We simulate precision measurements of gauge-mediated supersymmetry (SUSY) 
breaking (GMSB) parameters at a 500 GeV e+e" linear collider (LC) in the sce- 
nario where a neutralino is the next-to-lightest supersymmetric particle (NLSP). 
Information on the SUSY breaking and the messenger sectors of the theory is 
extracted from the measured sparticle mass spectrum and neutralino lifetime. 



GMSB0 is an attractive possibility for physics beyond the standard model (SM) 
and provides natural suppression of the SUSY contributions to flavour-changing 
neutral currents at low energies. In GMSB models, the gravitino G is the LSP with 

mass given by ttiq = ~ 2.37 (^joo^yv) '^here Vf is the fundamental 

SUSY breaking scale. The GMSB phenomenology is characterised by decays of 
the NLSP to its SM partner and the G with a non-negligible or even macroscopic 
lifetime. In the simplest GMSB realizations, depending on the parameters Afmoss, 
iVmcss, A, tan/3, sign{fi) defining the model, the JNLSP can be either the lightest 
neutralino Ni or the light stau fi. For this study □, we generated several thousand 
GMSB models following the standard phenomenological approach B and focused on 
the neutralino NLSP scenario, for which we selected several representative points 
for simulation. Our aim was to explore the potential of a LC in determining the 
GMSB parameters. Firstly, for a sample model with light sparticles, we considered 
a measurement of the spectrum via threshold scanning and found that iVmcss and A 
(Afmoss and tan/3) can be determined to 0.1 (1)% after a 200 fb^^ run with c.o.m. 
energy between 200 and 500 GeV. Then, we investigated A^i lifetime measurements 
in the whole allowed ctj^^ range, performing detailed event simulation for a set of 

representative GMSB models. Indeed, since the iVi lifetime is related to ^/F by 
IGttVf^ 1 / TOjv, 



IOCS \ 100 TeV 



V 100 GeV/ ' ^ ' 



the GMSB framework provides an opportunity to extract information on the SUSY 
breaking sector of the theory from collider experiments that is not available, e.g., 
in supergravity-inspired models. 

For the models we considered, the neutralino lifetime ranges from microns to 
tens of metres. While the lower bound on ct-^_^ comes from requiring perturbativity 
up to the grand unification scale i, the upper bound is only valid if the G mass 
is restricted to be lighter than about 1 keV, as suggested by some cosmological 
arguments □. This circumstance is summarised in Fig. |l|. 

"Work supported also by DESY, Hamburg. 
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Figure 1: Scatter plot of the neutralino NLSP lifetime as a function of the messenger scale Mmoss 
(a) and mjY^^ (b). For each set of GMSB model input parameters, we plot the lower limit on 

^ , corresponding to %/F ~ V^mcss = V A Afmoss • We use only models that fulfil the limit 

^ 1 keV => 

mess 

^ ^ 2000 TeV suggested by simple cosmology. 
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For given A^i mass and lifetime, the residual theoretical uncertainty on deter- 
mining \/F is due to the factor of order unity B in Eq. (|l|), whose variation is quite 
limited in GMSB models (cfr. Fig. 

For our simulations, and in particular for short cr^y^ measurements, it was 
fundamental to take the A^i Gf f decays into account, in addition to the dominant 
one Ni — > 7G. We performed a complete analysis of these channels and found that 
in most cases of interest for our study the total width is given approximately by 



r(iVi ^ //G) = r(7Vi 



^^^^ 



N 



2 In^ - - +r(7Vi ^ ZG)BiZ ^ //) 

n(2) 

where the expressions for the widths of the 2-body iVi decays are well-knownB. In 
Fig. ||b, the branching ratio (BR) of the Ni •jG decay is compared to those of 
Ni ZG and Ni — > /I'^G (in the on-shell approximation) and those of the main 
A^i — > //G channels (inclucUiig virtual-photon exchange contributions only). 

Using CompHEP 3.3. 18U together with a home-made lagrangian including the 
relevant gravitino interaction vertices in a suitable approximation, we also studied 
the kinematical distributions of the A'^i ffG channels and implemented the 
numerical results in our GMSB event generator. 

For this purpose, we used SUSYGEN 2.2/03 0, modified to include the 3-body 
neutralino decays as discussed above. For each sample GMSB model, we considered 
in most cases a LC run at a c.o.m. energy such that the only SUSY production 
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Figure 2: (a) Scatter plot showing the relation between the neutralino NLSP lifetime and the 
fundamental scale of SUSY breaking \/F, i.e. the factor B in Eq.(|l|), as a function of the neutralino 
mass in GMSB models of interest for the LC (100 frij^^ ^ 250 GeV). Big grey dots in represent 
models with a light ij-selectron (102—150 GeV), small black dots are for the heavier selectron case 
(150-430 GeV). (b) Scatter plot for the BR's of various TVi decay channels as a function of the A^i 
mass. Dots in different grey scale (colours) refer to the decays TVi — > ■yG, Ni — > ZG (including off- 
shell effects), and to hadrons or e+e^ plus gravitino via virtual photon, as labelled. For reference, 
we also report results for the 2-body ffi — > h^G decay in the on-shell approximation, whose BR 

is always negligible. 

process open is NLSP pair production e+e" NiNi, followed by decays through 
all possible channels. For more challenging models where the light SUSY thresholds 
are close to each other, we simulated also events from i?-slepton pair production and 
used some selection cuts to isolate the iViTVi events, for which the ATi production 
energy is fixed by the beam energy (we also took into account initial-state radiation 
as well as beamstrahlung effects), allowing a cleaner cr^^ measurement. 

The primary vertex of the events was first smeared according to the assumed 
beamspot size of 5 nm in y, 500 nru-in x and 400 /im in z and then the events were 
passed through a fulJ GEANT 3. 21 la simulation of the detector as described in the 
ECFA/DESY CDREI. The tracking detector components essential to our analysis 
included a 5-layer vertex detector with a point precision of 3.5 /im in r(j) and z, a 
TPC possessing 118 padrows with point resolution of 160 /xm in r4> and 0.1 cm in 
z. In addition we assumed an electromagnetic calorimeter with energy resolution 
given by (10.3/v^+ 0.6)%, angular pointing resolution of 50/\/E mrad and timing 
resolution of 2/y/E ns. The dimensions of the whole calorimeter (electromagnetic 
and hadronic) were 172 cm < r < 210 cm and 280 cm < \z\ < 330 cm. 

A single neutralino produced with energy Epf_^ will decay before travelling a dis- 
tance A with a probability given by P(A) = 1 — exp(— A/L), where L = CTfi_^ {Pi) Ni 



3 



'0 23 jO 40 bO 60 70 80 90 1 DO 

r (imn) 



(b) 



Statistical 



Calorimeter Timing 

< > 



Calorimeter Pointing 

< > 



3d tracl^ing 



Projective Tracking 



10 -0 



L (em) 



Figure 3: (a) Reconstructed projective radial distances, r, of the A^i — > ffG decay vertex, for 
a LC run on a short-lifetime model with L = 10 /im together with a fit to an exponential plus 
constant, (b) Summary of the techniques used here for a cr^^ measurement to 10% or better. 



is the Ni "average" decay length and iPj)ff^ = i^%J''^%^ ~ 1)"^^^- 

For L less than a few cm, we used tracking for measuring the vertex of Ni — > 
Gff decays. For the very short case, L less than a few hundred fim, the beamspot 
size becomes important and a 3D procedure is not appropriate. Instead, the re- 
constructed vertex was projected onto the xy plane, where the beamspot size is 
very small, and we used the resulting distributions to measure the Ni lifetime. We 
studied several GMSB models with cr^^ in the allowed range and found that the 
intrinsic resolution of the method was approximately 10 fim. An example of the 
reconstructed 2D decay length distribution for a challenging model where the neu- 
tralino lifetime can be very short □ is shown in Fig. ||a for statistics corresponding 
to 200 fb~^ (r is the xy component of A). 

For 500 fim < L < 15 cm, we used 3D vertexing to determine the decay length 
distribution and hence the lifetime of the Ni. Vertices arising from A^i 7G and 
photon conversions in detector material were essentially eliminated using cuts on 
the invariant mass of the daughter pairs together with geometrical projection cuts 
involving the mass of the Ni and the topology of the daughter tracks. Methods 
of measuring the Ni mass using the endpoints of photon energies or threshold 
techniques, together with details of the projection cuts have been describeda. Using 
200 fb~^ of data, we concluded that a lifetime measurement with statistical error 
of approximately 4% could be made using this method. 

For L larger than a few cm, we used the A'^i — > "/G channel. The calorimeter 
was assumed to have pointing capability, using the shower shapes together with 
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appropriate use of pre-shower detectors. Assuming the pointing angular resolution 
mentioned above, we demonstrateda how a decay length measurement can be made. 
We concluded that for lifetimes ranging from approximately 5 cm to approximately 
2 m this method worked excellently, with statistical precisions ranging from a few % 
at the shorter end to about 6% at the upper end of the range. We also investigated 
the use of timing information to provide a lifetime measurement but found it to be 
of less use than calorimeter pointing. However, the use of timing in assigning purely 
photonic events to bunch crossings and for rejecting cosmic backgrounds should not 
be underestimated. 

For very long lifetimes, we employed a statistical technique where the ratio 
of the number of one photon events in the detector to the number of two photon 
events was determined as a function of cr^y^ , requiring the presence of at least one 
displaced 7. This allowed a largely model-independent measurement out to cr^y^ ~ 
few lO's m. The possibility of using the ratio^f the number of no-photon SUSY 
events to one photon events was also discusseda. The latter allows a greater length 
reach, but relies on model-dependent assumptions. 

In Fig. we summarise the techniques we have used as a function of L for a 
sample model. The criterion for indicating a method as successful is a measurement 
of L and the iVi lifetime to 10% or better. It can be seen that L can be well measured 
for lO's of /im <, L lO's of m, which is in most cases enough to cover the wide 
range allowed by theory and suggested by cosmology. 

With reference to Eq. (|^), we note that a 10% error in cr^y^ corresponds to a 3% 

error in a/F. This is of the same order of magnitude as the uncertainty on the factor 
B, which parameterises mainly the different possible Ni physical compositions in 
GMSB models (cfr. Fig. ^). We also checked explicitly that, in comparison, the 
contributing error from a neutralino mass measurement using threshold-scanning 
techniques or end-point methods is negligible □. 

Hence we conclude that, for the models considered and under conservative as- 
sumptions, a determination of a/F with a precision of approximately 5% is achiev- 
able at a LC simply by performing Ni lifetime and mass measurements in the 
context of GMSB with neutralino NLSP. Less model dependent and more precise 
results can be obtained by adding information on the Ni physical composition from 
other observables, such as Ni decay BR's, cross sections, distributions. 
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